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Abstract 


Overdamping  that  is  common  to  a  single-degree-of-fteedom  damped  linear  vibratory  system 
was  extended  to  multidegree-of-freedom  damped  linear  system,  viz,  a  drive  train  situated  in  a 
typical  gas  turbine.  Inequalities  involving  the  mass,  damping,  and  stiffness  parameters  were 
derived  to  form  a  system  with  a  free  response  that  was  overdamped  in  each  respective  node.  A 
general  method  to  be  employed  in  establishing  the  design  parameters  for  designing  systems  to 
be  overdamped  in  each  mode  has  been  identified  for  purposes  of  analysis,  and  the  mediod  was 
utilized  to  a  four-degree-  of-freedom  model  of  a  drive  train  in  a  gas  turbine  engine  with  a  new 
look  at  a  solution  methodology  for  overdamping  considerations.  This  technique,  or  method,  for 
eliminating  oscillations  in  n-degree-of-freedom  lumped  parameter  systems  by  increasing  the 
amount  of  viscous  damping  in  the  system  has  been  illustrated  by  using  actual  data. 
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1.  Introduction 


Viscous  damping  can  limit  the  oscillations  in  mechanical  systems  such  as  in  drive  trains  of 
modern  tanks.  If  we  look  at  a  linear  one-degree-of-freedom  spring,  mass,  and  dashpot  arrangement, 
the  selection  of  the  proper  values  of  mass,  stiffness  and  damping  constants  to  produce  an 
overdamped  or  critically  damped  system  can  be  effectuated  without  too  much  difficulty.  The 
solution  of  a  constant  coefficient  second-order  ordinaiy  differential  equation  shows  that  if  c  ^  2  mk 
(where  ni,  c,  and  k  are  the  mass,  damping,  and  stiffness  coefficients,  respectively)  then  the 
mechanical  system  will  not  oscillate.  The  purpose  of  this  treatise  is  to  show  similar  inequalities  for 
nonoscillation  of  multidegree-of-freedom  systems,  especially  inherent  in  the  accessory  gear  box, 
reduction  gear  box  and  output  shaft  of  a  typical  gas  turbine  engine.  This  type  of  analysis  can  prove 
to  be  invaluable  for  the  diagnostics  and  prognostics  of  engines  as  illustrated  by  Helfman,  Dumer  and 
Hanratty  (1995). 

The  systems  considered  here  are  those  that  can  be  modeled  by  the  matrix  differential  equation 

Mx(t)  +  Cx(t)  +  Kx(t)  =  0,  (1) 

where  x  (t)  is  an  n-dimensional  vector  of  displacements  and  M,  C,  and  K  are  nxn  symmetric  matrices 
containing  the  physical  parameters  of  mass,  damping,  and  stiffness  constants.  It  is  further  assumed 
that  M  and  K  are  positive  definite  matrices,  and  that  C  is  at  least  a  positive  semi-definite  matrix.  The 
procedure  presented  in  this  work  will  take  advantage  of  the  derived  matrix  conditions  (Inman  and 
Andry  1980),  to  generate  nonlinear  algebraic  inequalities  for  the  physical  parameters  of  the  system. 
If  tlie  parameters  can  be  chosen  to  satisfy  these  inequalities,  the  resulting  transient  energy  response 
will  be  overdamped  in  each  respective  mode.  The  inequalities  are  stated  explicitly  in  temis  of  the 
mass,  damping,  and  stiffness  constants  of  the  system. 

The  exact  relations  for  overdamping  will  be  derived  for  a  general  two-degree-of-freedom  system. 
As  soon  as  the  damping  criterion  becomes  satisfied,  the  results  will  be  used  to  calculate  the 
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eigenvalues  of  the  system  to  illustrate  that  the  damped  system  is,  in  fact,  overdamped  in  each  mode. 
The  design  of  a  specific  four-degree-of-freedom  model  of  a  drive  train  in  a  gas  turbine  engine 
asscociated  with  the  reduction  gear  drive  train  in  a  tank  will  he  given  to  illustrate  the  problems 
encountered  in  more  practical  design  situations.  The  generalization  to  n  degrees  of  freedom  will 
become  obvious  from  these  examples. 

2.  Conception 


Inman  and  Andiy  (1980)  tell  us  that  if,  in  addition  to  the  restrictions  previously  listed,  the 
matrices  M,  C,  and  K  are  such,  that  the  matrix 


M-1/2  '  (2) 


is  positive  definite,  tlren  all  of  the  eigenvalues  of  equation  (1)  will  be  negative  real  numbers,  and 
hence  each  mode  will  be  classified  as  overdamped.  Since  M  is  symmetric  and  positive  definite,  it 
possesses  a  unique  positive  definite  square  root,  M'^"  with  inverse  M  Using  the  transformation 

X  =  y. 


equation  (1)  is  reduced  to 


y  +  Ay  +  By  =  0,  (3) 

where  A  =  and  B  =  M"''*  The  condition  for  overdamping  in  each  mode  for 

equation  (2)  is  that  the  matrix  A-2B'^  must  be  positive  definite.  Since  the  square  root  of  a  matrix 
is,  in  general,  harder  to  compute  then  the  square  of  a  matrix,  it  is  tempting  to  use  the  matrix  A^-4B 
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in  design  work.  Fortunately,  it  has  been  shown  (Bellman  1 968)  that  if  A--4B  is  positive  definite, 
then  so  is  A-2B'^l  Thus,  requiring  the  matrix  A--4B  to  be  positive  definite  ensures  that  each  mode 
of  equation  (2)  will  be  a  decaying  nonoscillating  (overdamped)  function  of  time. 

If  it  is  desired  to  make  the  solution  of  equation  (3)  overdaniped  in  each  mode  for  arbitrar}^  initial 
conditions,  then  it  suffices  to  choose  the  physical  constants  nij,  Cj,  and  k-,  so  that  A'-4B  is  positive 
definite. 

3.  Methodology 


To  illustrate  the  previously  mentioned  ideas,  consider  the  two-mass  arrangement  in  Figure  1. 
The  appropriate  matrices  for  the  equations  of  motion  are 


ni, 

0  ■ 

,  c  = 

C,  +  Cl  -C, 

,  and 

K  = 

k,  +  ki  -ki 

0 

m2 

-Cl  Cl 

-ki  ki 

Figure  1.  Two-Degrcc-of-Frecdom  System. 


The  matrix  M''^'  is 


m-'/2  ^ 


0 

l/^n^ 


(4) 


Denoting  the  element  of  a  generic  matrix  A  by  Ay  and  forming  the  matrix  A^-4B  yields 


(A--4B)„ 


(c,  +  c,)^ 

2 

m, 


Ct  k,  + 

+  — ^  -  4  — - ^ 

nij  nij 


Cj  C2  C2 

m,  ym,  m2 


4k^ 

- r -  +  - — 

m^ym,  nij  j/m, 


(6) 


mjm2 


(7) 


It  is  desired  to  choose  nij,  Cj,  and  k,  so  that  the  matrix  A^-4B  is  positive  definite.  A  necessary  and 
sufficient  condition  for  a  matrix  D  to  be  positive  definite  is  for  each  of  its  leading  principal  minors 
to  be  positive.  In  particular,  a  real  2x2  matrix  D  is  positive  definite,  if  and  only  if 

D,,>0,  and  D,|D22- D,2D2|>0  . 


Applying  these  inequalities  to  (A*  -  4B)  yields 


(c,  +03)^ 


>  4 


k,  +k2 


m,  1112 


m. 


(8) 


4 


X 


(c,  +c^y 


m, 


nij  m2 


-  4 


k,  +k^ 


m, 


C-, 


m 


2 

2 


O')  Iv-j 

+  — ^  -  4  — 


> 


/  2  2  \  ^ 

4k,  -  ^ 

^  ym,m2  m, ^111,1112 


(9) 


If  the  parameters  m|,  C;,  and  kj  are  now  chosen  to  satisfy  equation  (9)  (along  with  the  physical 
constraints  that  nij,  C|,  and  kj  are  all  positive)  then  equation  (1 )  will  be  overdamped  in  each  mode  and 
will  not  oscillate  when  perturbed  from  equilibrium.  In  total,  the  six  parameters  must  satisfy  eight 
inequalities  (Brent  1973;  Byrne  and  Hall  1973). 


The  approach  taken  here  was  simply  to  fix  the  values  of  nij  and  k;  and  chose  values  of  Cj  to  satisfy 
equation  (3).  For  example,  the  values 

m,  =  1  m,  =  1 
c,  =  4  C2  “  ^ 
k,  =  1  k2  =  2 


satisfy  equation  (9).  In  order  to  verify  that  this  set  of  values  implies  overdamping,  we  solve  the 
eigenvalue  problem  using  these  parameters.  This  yields  the  characteristic  polynomial 


A'*  +  14A-^  +  25A2  +  13A  +2  =  0 


with  eigenvalues 


X^  =-0.2891 , 
A2  =-0.4652, 


(10) 


A3  =-1.2389, 

and 

A4  =-12.0068. 

Thus,  the  design  procedure  yields  an  overdamped  response,  since  each  eigenvalue  is  a  negative  real 
number. 

4.  Design  Application 

The  process  described  here  may  be  useful  in  enhancing  the  survivability  of  certain  structures  by 
designing  them  to  have  an  overdamped  free  response.  In  order  to  illustrate  this  in  a  design  contest, 
we  consider  the  drive  train  of  a  gas  turbine  engine  (see  Figure  2).  The  numerical  values  for  inertia 
and  stiffness  are  listed  in  the  appendix,  along  with  the  definition  of  each  parameter. 


Figure  2.  Schematic  of  Turbine  and  Drive  Train  Components. 

In  order  to  produce  a  C  matrix  that  would  allow  inequalities  similar  to  equation  (9)  to  be 
formulated,  some  mechanism  must  be  available  for  adding  damping  to  the  system.  For  non- 
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rotational  s5^stems,  this  may  be  accomplished  by  the  use  of  shock  absorbers  or  linear  actuators.  For 
dampers,  this  may  be  useful.  Figure  2  indicates  the  addition  of  such  dampers  to  an  existing  system 
(i.e.,  c„  c,,  C3,  C5,  and  c^). 


and 


.T,  0  0  0 

0  J2  0  0 

0  0  J3  0 

0  0  0  .1, 


0 

1 

0 

0 

0 

1 

C,+C2 

-‘^2 

0 

0 

-c, 

J 

0 

0 

-C3 

+  d 

k,  -k,  0  0 

-k[  k, +k2  -k2  0 

0  -k2 

0  0  -k4n,n2  k^ 


(11) 


(12) 


(13) 


where  J,-  =  the  various  values  of  inertia;  C;  =  added  damping  constants;  d4  =  the  damping  constant, 
due  to  a  possible  connection  strut;  k,  and  k,  =  shaft  stiffness  constants;  and  kj  and  k4  =  stiffness 
constants  associated  with  the  transmission  and  gear  system,  fhe  transmission  has  a  gear  ratio  n,,  and 
the  timing  gear  has  a  ratio  of  n,.  Forming  the  matrix  A'-4B  yields 


(A^-4B)„ 


n 

^2 


.l,.l2 


(14) 
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(A2  -  4B),2  =  (A^  -  48).,  =  - 


C|(Ci  +c^)  ^  4k, 


j,yT^  J2\/X[I 


2  yhh 


(A2  -  4B),3  =  (A2  -  4B)3, 


Co 


(A^-4B),,  =  (A-4B)„  =  0, 


c,"  (c,  +  c,)-  c-T 

(A^-4B),,  =  -i-  *  ~  +  -  4 


J.J,  J,-  >A  ■'2 


■^2  \p2h  ■*:.\/'*2  '*3 


=  (A2-4B)3j, 


.  c. 

2  /m\  -  r  A  2  .iT^\  _ 


(A2-4B),^  =  (A"-4B> 


42 


C2  (c,  +  c.)^  Cj"  kj  ^  kj 
(A^-4B)33  =  - i-  -  4^ - 

‘'r's  T- 


J,-^  -'.vU 


(A^-4B)33  = 


43(03  *  C3)  _  03(0.  t  d,)  ^  ^ 

J3  ■'a  \/  '(l  -^-l 


=  (a2-4B),3, 


and 


(15) 


(16) 


(17) 


(18) 


(19) 

(20) 

(21) 


(22) 
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The  addition  of  c„  C2,  C3,  C5,  and  Cg  is  necessaiy  to  make  A--4B  positive  definite.  Requiring  the  four 
leading  principal  minors  of  the  4x4  matrix  A--4B  to  be  positive  yields  four  inequalities  in  the 
inertia,  damping  and  stiffness  parameters.  Using  the  values  lor  .T|  and  k|  (listed  in  the  appendix)  and 
choosing  q  to  satisfy  the  inequalities  yields 


c,  =  12.0000  X  10' 

N-m 

s/rad, 

c,  =  5.8653  X  10' 

N-m 

s/rad, 

C3  =  1.4700  X  10- 

N-m 

s/rad, 

c,  =  3.5300  X  10- 

N-m 

s/rad, 

C5  =  18.0000  X  10' 

N-m 

s/rad, 

Cg  =  7.1347  X  10' 

N-m 

s/rad 

as  one  possible  solution  for  the  added  damping  constants. 

The  characteristic  polynomial  for  this  system  is 

2.024397  X  10'^  11*  + 8.57279  x  10'^  V+  1.186625 
+  58.521594  +  7.830627  x  10- 

+  2.5788556  X  10^  A.^  +  5.7545081  x  10^  X’ 

+  5.7354329  x  lOU  -i- 16.620975  =  0,  (24) 

which  has  roots 
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A,  =  -0.002, 

X,  =  -1.707, 

A3- -14.814, 

A,  =-68.894, 

Aj  =-1.5  X  10% 

Ag  =  -1.85  X  10% 

A7  = -0.943  +2.2161, 

and 

As  =  -0.943  -2.2161 


5.  Conclusion 


A  method  accompanied  by  the  complexity  of  the  process  and  its  level  of  applicability  has  been 
presented  for  eliminating  oscillations  in  n-degree-of-freedom  lumped  parameter  systems  by 
increasing  the  amount  of  viscous  damping  in  a  four-degree-of-freedom  system. 

Another  method  available  to  produce  total  overdamping  is  given  in  Beskos  and  Boley  (1980)  for 
two-degree-of-freedom  systems.  Unfortunately,  to  extend  the  process  in  Beskos  and  Boley  (1980) 
to  n  degrees  of  freedom  requires  a  closed-form  solution  of  polynomials  of  degree  (n  -  1).  However, 
the  method  here  requires  only  the  numerical  solution  of  nonlinear  equalities.  In  addition,  tor  the 
two-degree-of-freedom  case,  the  method  presented  in  Bellman  (1970)  allows  only  the  parameters 
c,  and  Cj  to  be  adjusted.  As  an  alternative,  the  method  presented  here  allows  all  of  the  parameters 
mj,  Cj,  and  kj  to  be  sensitized  i.e.,  adjusted  to  effectuate  mass,  damping  and  stiffness  characteristics. 
Hence,  it  seems  quite  apparent  that  this  method  is  more  advantageous  for  design  work  especially 
when  it  involves  computationally  intensive  operations. 
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Appendix: 
Parameter  Values 
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The  parameter  values  of  the  turbine  for  the  case  of  tip  speed  (ratio  of  2)  and  turbine  rotational 
speed  of  5,000  RPM  are 

J,  =  J,  =  1/2  of  turbine  motor  inertia  =  1 .65  x  10“*  N  -  s''  -  m  (1 .46  x  lOMb  -  s^  -  in) 

J3  =  transmission  inertia  =  2.43  x  10‘  N  -  s"  -  in  (2.1 5x10'’  lb  -  s'  -  in) 

J4  =  generator  inertia  =  3.06  N  -  s"  -  in  (27.1  lb  -  s"  -  in) 

n,  =  transmission  gear  ratio  =  35.6 

„  178,000  ^000 

=  pulley  gear  ratio  =  =  0.999 

kj  =  rotor  stiffness  =  1.65  x  10^  N-m/rad  (1.46  x  10^  lb  ~  iii/rad) 

k2  =  shaft  stiffness  =  2.69  x  10^  N  -  ni/rad  (2.39  x  10^  lb  -  in/rad) 

Kj  =  transmission  shaft  stiffness  =  1.41  x  10^  N  -  m/rad  (1.25  x  10^  lb  -  in/rad) 

Ki.|  =  generator  shaft  stiffness  =  2.1 0  x  10'*  N  -  m/rad  (1 .86  x  1 0*^  lb  -  in/rad) 
n  ^  n  ^  K  K 

=  _J — I — ! — H  =  2.62  X  10^  N  -  m/rad  (2.39  x  10'^  lb  -  in/rad) 
k,  +  n2  K„ 

^  =  2.07  X  10^11  -  m/rad  (1 .83  x  10''’  lb  -  in/rad) 

11,-112 
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